From here:
postgibbspostgibbsgibbs_samples <-
purrr::map2_dfr(.x = rep(c(1:10),
times = 7),
.y = rep(c("1", "2", "3alt", "5", "7", "8", "9"),
times = 10),
~ read_gibbs_samples(iteration = .x,
region = .y)) %>%
select(iter, dataset, everything(), -X1, -X3) %>%
purrr::set_names("iter", "dataset", "round", c(1:14)) %>%
tidyr::pivot_longer(cols = c(`1`:`14`),
names_to = "param") %>%
mutate(param = case_when(param == "1" ~ "dir1dir1",
param == "2" ~ "dir1dir2",
param == "3" ~ "dir1mat1",
param == "4" ~ "dir1mat2",
param == "5" ~ "dir2dir2",
param == "6" ~ "dir2mat1",
param == "7" ~ "dir2mat2",
param == "8" ~ "mat1mat1",
param == "9" ~ "mat1mat2",
param == "10" ~ "mat2mat2",
param == "11" ~ "mpe1mpe1",
param == "12" ~ "mpe2mpe2",
param == "13" ~ "res1res1",
param == "14" ~ "res2res2"),
iter = as.character(iter))
gibbs_mce <-
purrr::map2_dfr(.x = rep(c(1:10),
times = 7),
.y = rep(c("1", "2", "3alt", "5", "7", "8", "9"),
times = 10),
~ read_gibbs_mce(iteration = .x,
region = .y)) %>%
left_join(gibbs_samples %>%
mutate(iter = as.integer(iter)) %>%
group_by(iter, dataset, param) %>%
tally(name = "n_samples")) %>%
select(iter, dataset, param, n_samples, everything())
gibbs_psd <-
purrr::map2_dfr(.x = rep(c(1:10),
times = 7),
.y = rep(c("1", "2", "3alt", "5", "7", "8", "9"),
times = 10),
~ read_gibbs_psd(iteration = .x,
region = .y)) %>%
left_join(gibbs_samples %>%
mutate(iter = as.integer(iter)) %>%
group_by(iter, dataset, param) %>%
tally(name = "n_samples")) %>%
select(iter, dataset, param, n_samples, everything())
plot_gibbs_iter <-
function(df) {
df %>%
filter(!stringr::str_detect(param, "res|mpe")) %>%
ggplot(aes(x = round,
y = value,
color = param)) +
geom_line() +
ggsci::scale_color_ucscgb() +
theme_classic() +
facet_wrap(~ dataset, nrow = 3)
}
plot_gibbs_region <-
function(df) {
df %>%
filter(!stringr::str_detect(param, "res|mpe")) %>%
mutate(iter = as.numeric(iter)) %>%
arrange(iter) %>%
mutate(iter = glue("Iteration {iter}"),
iter = forcats::fct_inorder(as.factor(iter))) %>%
ggplot(aes(x = round,
y = value,
color = param)) +
geom_line() +
theme_classic() +
ggsci::scale_color_ucscgb() +
labs(x = "Round",
y = "Value",
color = "Parameter") +
facet_wrap(~ iter, ncol = 2)
}
gibbs_samples %>%
filter(dataset == "3v1") %>%
plot_gibbs_region()
gibbs_samples %>%
filter(dataset == "3v2") %>%
plot_gibbs_region()
gibbs_samples %>%
filter(dataset == "3v3alt") %>%
plot_gibbs_region()
gibbs_samples %>%
filter(dataset == "3v5") %>%
plot_gibbs_region()
gibbs_samples %>%
filter(dataset == "3v7") %>%
plot_gibbs_region()
gibbs_samples %>%
filter(dataset == "3v8") %>%
plot_gibbs_region()
gibbs_samples %>%
filter(dataset == "3v9") %>%
plot_gibbs_region()